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ABSTRACT 

We have examined whether the major axes of rich galaxy clusters tend to point toward their nearest 
neighboring cluster. We have used the data of Ulmer, McMillan, and Kowalski, who used position 
angles based on X-ray morphology. We also studied a subset of this sample with updated positions 
and distances from the MX Northern Abell Cluster Survey (for rich clusters (R > 1) with well known 
redshifts). A Kolmogorov-Smirnov (KS) test showed no significant signal for nonrandom angles on any 
scale < 100/i _1 Mpc. However, refining the null hypothesis with the Wilcoxon rank-sum test, we found 
a high confidence signal for alignment. Confidence levels increase to a high of 99.997% as only near 
neighbors which are very close are considered. We conclude there is a strong alignment signal in the 
data, consistent with gravitational instability acting on Gaussian perturbations. 

Subject headings: galaxies: clusters: general — large-scale structure of universe 
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1. INTRODUCTION 

It is well documented that clusters of galaxies tend to 
be elongated, elliptical systems (eg., Carter & Metcalfe 
1980) giving them major axes and "position angles" in 
the sky. Bingelli (1982) found that the major axes of 
rich galaxy clusters have the tendency (in projection) to 
point toward their nearest neighbor cluster whose dis- 
tance, d n , was closer than ~ 30h~ 1 Mpc. Since then, 
there have been multiple studies on whether this "Bingelli 
effect" actually exists. Much of this literature supports 
the reality of the effect. Flin (1987) and Rhee & Kat- 
gert (1987) both found significant alignments for d n less 
than ~ 30/i _1 Mpc. West (1989) and Rhee, van Haarlen 
& Katgert (1992) have also detected a signal for align- 
ment. Plionis (1994) found weak alignment signals up to 
d n r~j 60/i -1 Mpc, with more significant cluster alignments 
on smaller scales ( 10 — 30/i~ 1 Mpc). Moreover, West, Jones 
and Forman (1995) found evidence that galaxy cluster sub- 
structure tends to be aligned with its host cluster and sur- 
rounding environment out to ~ lOft-^Mpc, which might 
help explain these alignments. 

Not all authors favor an alignment effect, however. Both 
Struble & Peebles (1985) and Ulmer, McMillan & Kowal- 
ski (1989) (hereafter UMK) found no significant evidence 
that clusters point toward their nearest neighbor (see how- 
ever, Argyles et al. 1986). 

Galactic positions may not be good tracers of the shape 
of a cluster, for galaxies contribute discreteness noise. 
Most clusters contain much more mass in hot, X-ray emit- 
ting gas than that of the galaxies themselves. Dark matter 
contributes more mass to the system than gas and galaxies 
combined. Thus, the shape of the actual cluster mass can- 
not be directly seen. However, it is believed that the X-ray 
emitting gas within a cluster traces its gravitational poten- 
tial (Sarazin 1986). X-ray morphology is, then, probably 
the best observable for determining galaxy cluster shape 
and orientation. 

Cluster alignments are not crucial in distinguishing cos- 



mological models, but they are additional evidence in sup- 
port of the gravitational instability hypothesis of structure 
formation (Shandarin & Klypin 1984; Splinter et al. 1997; 
Onuora & Thomas 2000). 

For these reasons, the negative results of UMK are in- 
teresting. Whereas most authors used galaxies to define 
ellipticity and spatial orientation, both UMK and West et 
al. chose to use X-ray morphology. UMK did not find 
a statistical alignment for any nearest neighbor distance 
scale. There are many more papers (including West et 
al.) that found an alignment effect than didn't. Since 
UMK used the shape of the X-ray gas in their search 
for alignment, their negative results are even more impor- 
tant. Most analyses of numerical simulations of structure 
formation by gravitational instability from Gaussian ini- 
tial perturbations predict alignments on some scale, which 
provides some physical motivation for detecting the align- 
ments searched for by UMK. 

Both Onuora & Thomas 2000 and Splinter et al. 1997, 
for example, predicted alignments for d n of at least 
15/i _1 Mpc for standard CDM models. Onuora & Thomas 
extended this distance to 30ft _1 Mpc for ACDM models. 
Both these studies showed that the predicted alignment 
differences for individual cosmological background models 
are not practical means for determining cosmological pa- 
rameters, such as ft. However, Splinter et al. presented 
evidence that cluster ellipticity and the scale dependence 
of cluster alignments probe the primordial power spectrum 
independent of the parameters of the background cosmol- 
ogy. The alignments in these simulations fit a general 
picture of cluster formation by hierarchical clustering in 
which material falls into the cluster along the large scale fil- 
amentary structure, as interpreted by Shandarin & Klypin 
(1984). This picture has been supported by dynamical ev- 
idence of drainage along such filaments (Novikov et al. 
1999). 
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2.1. Subject Clusters and Position Angles 

UMK determined the major axis orientation of 46 X-ray 
clusters observed by Einstein (UMK Table 1). The ma- 
jor axis served to define a position angle on the celestial 
sphere (measured counter-clock- wise from north). UMK 
used both R = and R > 1 Abell clusters in their origi- 
nal analysis. However, R = clusters were never part of 
Abell's (1958) statistical sample and should not be used 
in nearest-neighbor analyses (since one needs to be sure 
of the existence of both the source cluster and its nearest 
neighbor) . UMK created smaller subsets of clusters (out of 
their original 46) for analysis as well. We will re-examine 
only the largest (46 cluster) subset from UMK, since this 
was the only one which was large enough to use after we 
eliminated clusters for which we did not have good red- 
shift information. We analyzed the data in UMK Table 
1 as given, again with updated redshifts (results we show 
here), and then re-analyzed a subset of it which met our 
stringent selection criteria as described below. 

2.2. Potential Neighbor Sample 

We start out with the 46 clusters from UMK Table 1. 
However, over the past ten years many clusters have new 
or revised redshifts, and so we also searched an updated 
Abell cluster redshift survey to revise any of the nearest 
neighbors in UMK and also to apply constraints that will 
limit any selection effects and biases in our data. Our 
angles are measured counter-clockwise from North as in 
UMK. Our distances are measured for a Friedmann Uni- 
verse with = and h = i?o/100km s _1 Mpc _1 . 

The recent success of galaxy cluster surveys (Miller et 
al. 1999; 2000 and references therein) gives us a large and 
uniform sample, as compared with previous efforts. The 
cluster redshifts and coordinates we have used are mainly 
from the MX Northern Abell Cluster Survey (MX) (Slin- 
glend et al. 1998; Miller et al. 2000). There are 256 Abell 
R > 1 clusters (i.e. dec. > —27°) having measured red- 
shifts and within 0.012 < z < 0.10 and \b\ > 30°. These 
clusters have an average of 25 measured galaxies each, and 
90% have more than one measured redshift. Having mul- 
tiple redshifts per cluster is quite important, as the elon- 
gation in redshift space due to internal velocity disper- 
sion in the cluster is easily comparable to typical nearest 
neighbor distances. Miller et al. (1999) found that cluster 
redshifts based on only one galaxy redshift are erroneous 
by ±500km s _1 41% of the time. While this error may 
seem small (<~ 5/i _1 Mpc), such an error can easily throw 
off the determination of the nearest-neighbor in dense en- 
vironments. This northern hemisphere sample is nearly 
complete to z = 0.10 (e.g. Miller & Batuski 2000). To 
define potential nearest neighbors, we searched through 
these 256 clusters noting the distance and direction to the 
nearest neighbor and the distance to the nearest edge of 
the survey. We then made a potentially important cut: 
If the boundary of the survey region were closer than the 
nearest neighbor, the pair would not be used for our align- 
ment statistics. This is because there is a potential nearer 
neighbor hidden outside the boundary. 

After we applied the above constraints to the original 
UMK data, we were left with 25 clusters out of their origi- 
nal 46. We call this our Statistical Sample and we present 
our data in Table 1. We also found that ten of these 25 



clusters now have different nearest-neighbors than those 
found by UMK. In four of these ten cases, UMK used an 
R = cluster as their nearest neighbor. In other words, 
the original UMK dataset remains relatively unchanged, 
with 40/46 clusters having the same nearest-neighbor in 
our new analysis. It is worth noting that the constraints 
that we apply to the original UMK data prevent a large 
number of biases from entering the analysis. For example, 
by using only R > 1 clusters we are ensuring that our base 
cluster subset is a statistically complete sample. We also 
use clusters with multiple-galaxy determined redshifts so 
that we can be sure of their location. Finally, perhaps our 
most important improvement over the UMK data, is that 
we exclude clusters where the edge of the survey is closer 
than the nearest-neighbor. 

2.3. Analysis 

We first analyze the UMK data as given (in their Ta- 
ble 1) with revised nearest neighbors for six clusters. Our 
purpose in re-analyzing the original UMK data is to de- 
termine whether their statistical analyses were sensitive 
enough to actually detect a cluster alignment. We note 
that the UMK dataset contains R = clusters and does 
not exclude clusters that have a survey-edge which is closer 
than the nearest neighbor. We then analyze our subset of 
25 clusters which meet our stringent selection criteria de- 
signed to produce a uniform well-controlled sample. 

UMK have provided a position angle for the major axis 
of the cluster as measured via the X-ray emission. As this 
is an orientation, not a direction, the angle between it and 
the projected direction to the nearest neighbor, the point- 
ing angle <f> p can only have a range —90 < <j) p < 90. We 
also assumed that the sign is not significant, so we examine 
< \4> p \ <90. This coincides with the UMK procedure. 

We therefore define alignment as a tendency for the an- 
gles \4>p\ to be smaller than they would be if distributed 
isotropically, that is, uniformly over this interval. Most 
previous work has not really tested for alignment - rather, 
it has tested for any kind of anisotropy. 

2.3.1. Re-analysis of UMK Data 

We repeated the Kolmogorov-Smirnov (KS) test as used 
by UMK to test against a distribution sampled from a 
population uniformly distributed over this interval, as <f> p 
would be if there are no correlations. The KS test (Lehman 
1975) is used to test against the null hypothesis that the 
sample (our angles cj> p ) could be drawn from a parent pop- 
ulation of 4> p . We have no a priori reason to believe that 
the parent distribution of p 's should be anything other 
than random. In other words, the null hypothesis for the 
KS test is that our observed sample is the same as our 
parent (random) distribution. The KS test measures the 
significance of the null hypothesis being false. Although 
widely used in the "cluster alignment" literature, it in fact 
is a test for non-uniformity (since our hypothetical parent 
population has a uniform distribution of <fi p ). The KS test 
uses the maximum value of the difference between cumu- 
lative distribution functions as its diagnostic. Thus, KS is 
excellent at detecting any deviation of the sample from the 
parent population. However, this very attribute makes it 
weaker at detecting a more specific signal (in our case the 
signal of alignment) than some other tests. Using the KS 
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test, our results agree with UMK. We find that the null hy- 
pothesis can be discounted with only 73% confidence. This 
is too small a confidence level to rule out the null hypoth- 
esis. We have also looked only at those clusters with d n 
below some critical value, d c . For d c down to 10ft- _1 Mpc, 
the confidence for ruling out the null hypothesis is still too 
small. 

Refinement of the null hypothesis makes an enormous 
difference. We are not looking for just any difference, we 
are looking for alignments, which means the angles <f) p are 
systematically lower than they would be if drawn from 
a uniform parent population of <f>. The Wilcoxon rank- 
sum test (WRS) (Lehmann f 975) tests for this. The null 
hypothesis of WRS is that the sample is not systemati- 
cally smaller or larger than the parent population. Thus, 
while WRS is more sensitive to alignment (small angles), 
it would not (for example) be sensitive to a tendency for 
the pointing angles to clump around 45 degrees. The WRS 
ranks the populations and is thus sensitive to alignment 
differences between the real and assumed parent popula- 
tion. The WRS returns a signed result, indicating that 
our sample has systematically higher or lower (ftp's. By 
defining more precisely what is tested, a great increase in 
statistical power can be achieved. 

Using WRS on the updated UMK data, we found no 
significant effect. However, if we use it only on the UMK 
clusters with d n < 30/i _1 Mpc, we find (for the 38 UMK 
clusters that meet this criterion) 99.8% confidence in align- 
ment. Restriction to smaller distances produces align- 
ments with much higher confidence, as seen in Table 2. 
Refinement of the null hypothesis combined with restric- 
tions to nearer neighbors produces a strong signal. We also 
note that the original UMK data (i.e. with six erroneous 
nearest neighbors) also shows a similarly strong signifi- 
cance (99.6%) for alignment. In other words, had UMK 
used a more targeted statistical tool (such as the WRS), 
they too would have certainly noted this alignment effect. 

2.3.2. Analysis of Revised UMK Data 

We next examined the set of our 25 clusters with cor- 
rected distances in our Statistical Sample. We observe 
that this set has fractionally fewer nearest neighbors at 
large distances than in UMK's Table 1, suggesting correc- 
tions of misidentification of near neighbors in that. How- 
ever, we still found no significant signal with the KS test. 
Applying WRS to our entire controlled sample of 25 clus- 
ters produced, as it did with the entire UMK sample, no 
substantial confidence in alignment. Restricting the study 
to the clusters with d n < 20/i _1 Mpc leaves 17 clusters, 
with a confidence level of 98.4%. More restrictive limits 
on distance produces increased confidence, up to very high 
levels, as with the UMK data. 

Table 2 shows the results of the WRS test applied to 



both the original UMK data (with six revised cluster 
neighbors) as well as to the 25 clusters with corrected dis- 
tances and stringent constraints to account for any bias- 
ing or selection effects. Table 2 shows the strength of the 
alignment effect as we vary the restrictions on how near 
the neighbors must be in order to be considered. 

There are suggestive trends in the data. Although the 
UMK data contain nearly twice as many clusters, as more 
restrictive distance cuts are applied, the number of surviv- 
ing clusters converge until both sets are nearly the same 
size for d n < 10/i _1 Mpc. Because clusters have correlated 
spatial positions, a random error in cluster position mea- 
surement is much more likely to move a cluster away from 
a near neighbor than to create a spurious one. Our more 
well-controlled sample has a smaller mean nearest neigh- 
bor distance. 

The UMK result, in this case, nevertheless shows a 
strong signal when even moderately restrictive cuts are 
made. By making these cuts we are most likely remov- 
ing erroneous near neighbor identifications. The surviving 
close near neighbors in UMK may well be correct, in spite 
of being at low galactic latitude (where obscuration is a 
problem), being in poorly sampled regions, or being too 
close to a survey boundary. 

3. CONCLUSIONS AND DISCUSSION 

Re-examining the UMK Einstein cluster data, we find 
that testing for alignment (with WRS) rather than for 
any departure from uniformity of angles (as with KS) 
allows us to find a significant signal for alignment in 
this data for clusters with nearest neighbors at distances 
d n < 30/i _1 Mpc. When we use a more stringently defined 
sample, we still find a strong signal for alignment, reaching 
3.34(7 when we restrict d n < 10/i~ 1 Mpc. Thus, refinement 
of the null hypothesis has proven crucial in finding the 
alignment signal. Use of the well-controlled sample was 
in this case not necessary to find a signal, but it confirms 
(with a lower confidence due to a smaller sample size) that 
the X-ray emission from galaxy clusters does tend to point 
to the nearest cluster neighbor. 

A potential weakness in most alignment studies is the 
search for nearest neighbor alignment rather than super- 
cluster axis alignment. In the future, we plan to examine a 
larger sample using the supercluster axis finding procedure 
defined in Novikov et al. (1999). The work we present here 
serves to remove an inconsistency from the data analysis 
and help clarify some issues of statistical methodology. 

ALM wishes to acknowledge support from the National 
Science Foundation under grant number AST-0070702, the 
University of Kansas General Research Fund and the Na- 
tional Center for Supercomputing Applications. 
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Table 1 

DATA SUMMARY FOR STATISTICAL SAMPLE 







Cluster 










Nearest 


Neighbor 






Abell no. 


R.A. 


dec. 


Dist. 


<t> 


Abell no. 


R.A. 


dec. 


Dist. 


dn 


4>v 


85 


0.652 


-9.617 


162 


166 


87 a 


0.675 


- 10.067 


161 


2.1 


23.6 


119 


0.897 


-1.533 


130 


18 


168 


1.210 


-0.017 


132 


11.5 


54.1 


154 


1.138 


17.400 


185 


57 


150 


1.110 


12.900 


172 


19.6 


51.7 


168 


1.210 


-0.017 


132 


156 


119 


0.897 


-1.533 


130 


11.5 


83.9 


399 


2.920 


12.817 


210 


31 


401 


2.937 


13.383 


214 


4.3 


7.8 


401 


2.937 


13.383 


214 


38 


399 


2.920 


12.817 


210 


4.3 


13.8 


1367 .... 


11.698 


20.117 


65 


137 


1656 .... 


12.957 


28.250 


69 


22.3 


70.3 


1656 .... 


12.957 


28.250 


69 


49 


1367 .... 


22.330 


11.698 


65 


22.3 


17.7 


1767 .... 


13.570 


59.467 


204 


145 


1904 .... 


14.338 


48.783 


205 


44.9 


12.2 


1775 


13.660 


26.617 


208 


117 


1795 b 


13.778 


26.833 


184 


24.8 


34.0 


1795 .... 


13.778 


26.833 


184 


21 


1831 b . . . . 


13.948 


28.233 


179 


9.7 


40.2 


1809 .... 


13.847 


5.400 


228 


1 


1780 a .... 


13.702 


3.133 


227 


12.4 


42.8 


1904 .... 


14.338 


48.783 


205 


84 


1767 .... 


13.570 


59.467 


204 


44.9 


48.8 


1991 .... 


14.870 


18.833 


171 


56 


1913 .... 


14.408 


16.900 


154 


25.7 


18.4 


2029 .... 


15.142 


5.950 


224 


131 


2028 .... 


15.118 


7.717 


225 


7.1 


37.5 


2063 .... 


15.343 


8.817 


104 


55 


2147 .... 


16.000 


16.033 


103 


21.7 


1.2 


2065 .... 


15.343 


27.900 


210 


151 


2056 a . . . . 


15.285 


28.450 


216 


7.0 


28.7 


2079 .... 


15.433 


29.050 


192 


41 


2092 .... 


15.522 


31.317 


194 


8.9 


10.5 


2107 .... 


15.627 


21.933 


121 


167 


2152 b . . . . 


16.052 


16.583 


121 


17.0 


37.0 


2124 .... 


15.718 


36.217 


192 


135 


2122 a . . . . 


15.710 


36.283 


192 


2.7 


16.2 


2147 .... 


16.000 


16.033 


103 


159 


2151 a . . . . 


16.050 


17.883 


108 


5.9 


43.1 


2151 .... 


16.050 


17.883 


108 


116 


2147 a . . . . 


16.000 


16.033 


103 


5.9 


86.1 


2152 .... 


16.052 


16.583 


121 


113 


2151 .... 


16.050 


17.883 


108 


13.0 


65.7 


2199 .... 


16.448 


39.633 


88 


43 


2197 .... 


16.442 


41.017 


91 


3.4 


46.7 


2670 .... 


23.860 


-10.683 


221 


104 


2659 a . . . . 


23.708 


-15.750 


226 


22.0 


79.8 



Note. — Key to columns is as follows: 
Col. (1). — Abell cluster number. 

Cols. (2)-(3). — Right ascension (hours) and declination (degrees) of the cluster (1950 epoch). 

Col. (4). — Distance to the cluster (h~ Mpc). 

Col. (5). — Position angle of the cluster. 

Cols. (6) - (10) — For the nearest neighbor cluster. 

Col. (6). — Number of the Abell cluster. 

Cols. (7) - (8). - Right ascension (hours) and declination (degrees) of the cluster (1950 epoch). 
Col. (9). — Distance to the cluster (h _1 Mpc). 

Col. (10). - Distance between the cluster (1) and its nearest neighbor (6) (A _1 Mpc). 
Col. (11). - "Pointing" angle in degrees; see text for definition. 
a Different nearest neighbor than that found by UMK. 

6 Different nearest neighbor that that found by UMK. However, UMK used an R = cluster. 



Table 2 
Wilcoxon Rank-Sum Results 



(fr^Mpc) 


UMK Data (rev.) 


Stat. Sample 




Sig. 




Sig. 




10 


4.00(7 (99.997%) 


12 


3.34c- (99.96%) 


11 


20 


3.45er (99.97%) 


29 


2.14(7 (98.4%) 


17 


30 


2.90<7 (99.8%) 


38 


1.27(7 (89.8%) 


23 


None 


0.94(7 (82.6%) 


46 


0.87(7 (80.8%) 


25 



5 




20 40 60 

Position Angle 

Fig. 1. — The distribution of position angle against nearest neighbor distance is shown for the UMK original 46 clusters. Six of these 
clusters have revised nearest neighbor distances and pointing angles. 




40 

Position Angle 

Fig. 2. — The distribution of position angle against nearest neighbor distance is shown for our sample of 25 clusters, using UMK position 
angles with corrected distances than were found in UMK. 
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